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Abstract 

We investigate the radiative effects of dust aerosols in the NASA GEOS-5 atmospheric 
general circulation model. A sectional aerosol and cloud microphysics model (CARMA) is 
included in GEOS-5. CARMA treats the dust aerosol lifecycle, including sources, transport 
evolution, and sinks. The CARMA aerosols are radiatively coupled to GEOS-5, and we 
perform a series of AMIP-style simulations in which dust optical properties (spectral 
refractive index and particle shape) are varied. Spherical dust optical properties are from 
Mie theory, while non-spherical properties are drawn from a recent database for tri-axial 
ellipsoids. The simulated dust distributions generally compare well to data from the space- 
based MODIS, MISR, and CALIOP, the ground-based AERONET, and surface measurements 
of dust deposition and mass concentration. We show significant variability in simulated 
summertime Saharan dust distributions resulting from different choices of dust optical 
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properties. Dust atmospheric heating enhances surface winds over important dust sources, 
increasing emissions. We find that increased dust absorption leads to a strengthening of 
the summertime Hadley cell circulation, increasing dust lofting to higher altitudes, and 
strengthening the African Easterly Jet. This leads to a longer atmospheric residence time 
and generally more northward transport of dust in simulations with the most absorbing 
dust optical properties, in best agreement with observations. We find that particle shape 
has a minor effect on dust transport and distribution, although total atmospheric forcing is 
enhanced by as much as 30% over land for simulations incorporating a spheroidal shape 
distribution versus ellipsoidal or spherical shapes. 

1. Introduction 

Mineral dust aerosols play many roles in Earth’s climate system. Scattering and absorption 
of shortwave (i.e., solar) and longwave (i.e., thermal) radiation by dust particles directly 
influence Earth’s radiation balance (Sokolik and Toon, 1996; Balkanski et al., 2007). The 
indirect effect of dust particles on climate includes their modification of cloud properties. 
For example, by acting as an additional source of cloud condensation nuclei (CCN), dust 
particles can lead to smaller cloud droplets, hence limiting droplet coalescence efficiency 
and suppressing precipitation (Rosenfeld et al., 2001). The facility of dust particles to act as 
ice nuclei in clouds offers another pathway toward indirect climate effects, where dust is 
observationally linked to glaciation in mildly super-cooled altocumulus clouds (Sassen et 
al., 2003). Dry air layers associated with Saharan dust outbreaks may play a role in 
suppressing convection and tropical cyclone formation in the Atlantic (Dunion and Velden, 
2004; Wong and Dessler, 2005), and satellite observations suggest an inverse relationship 
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between the presence of dust and tropical cyclone activity in the same region (Evan et al., 
2006; Lau and Kim, 2007). Observational evidence from meteorological analyses suggests 
a role for radiative heating of Saharan dust layers in modulating the amplitude of easterly 
waves in which the dust is typically transported (Jones et ah, 2004). Biological productivity 
(i.e., photosynthetic activity in phytoplankton) in many ocean regions is limited by the 
availability of iron, which has as a major and variable source the deposition of mineral dust 
aerosols (Martin, 1990; Fung et ah, 2000; Mahowald et ah, 2009). To the extent that an 
enhanced supply of nutrients from dust deposition mitigates the iron stress limitation in 
these ecosystems, dust can modulate the oceanic uptake of C0 2 at the ocean-air interface 
(see review in Jickells et ah, 2005). Similarly, dust can be an important nutrient source to 
terrestrial ecosystems (Swap et ah, 1992). Atmospheric dust aerosols play a role in 
tropospheric photochemistry pathways by altering photolysis rates and acting as a surface 
on which heterogeneous reactions can occur (Dentener et ah, 1996; Bian and Zender, 2003; 
Bauer et ah, 2004). Finally, long-range transport of dust aerosols are a source of fine 
particulate matter associated with poor air quality episodes (Prospero, 1999) and are a 
vector for transporting microorganisms associated with coral reef mortality (Shinn et ah, 
2000 ). 

Because of the importance of mineral dust aerosols in Earth’s climate system there has 
been a considerable effort to model their temporally and spatially varying distributions in 
chemical transport models and global climate models (Tegen and Fung, 1994; Mahowald et 
ah, 1999; Ginoux et ah, 2001; Zender et ah, 2003; Colarco et ah, 2003; Su and Toon, 2009; 
Nowottnick et ah, 2010). Synthesis studies of various modeling efforts reveal the 
importance of dust to the total aerosol loading, but also show wide variability in modeled 
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emissions and atmospheric burden (Zender et al., 2004; Textor et al., 2006; Huneeus et al., 
2011). There also remain large differences in the treatment of the particle size distribution 
and its evolution in dust models that can greatly affect the dust radiative impacts (Kok 
2011). These issues are convolved with model-dependent factors that impact source, sink, 
and transport terms in the aerosol evolution (e.g., spatial resolution, dynamical core, 
turbulent and convective parameterizations, precipitation, etc.). 

The direct radiative effects of dust aerosols on Earth’s climate thus remain uncertain in 
both magnitude and even sign in recent assessments (IPCC 2007). In addition to the 
uncertainties in atmospheric burden and particle size distribution noted above, there 
remains an uncertainty in the dust optical properties themselves. These are manifest in the 
composition of the dust — generally a function of source and particle size — and these will 
have important effects on the radiative properties of the aerosol (Sokolik and Toon, 1999; 
Quijano et al., 2000). The typical modeling paradigm is to carry the size resolved mass of 
dust as a number of tracers, and for each tracer to specify optical properties that are a 
function of size and the spectrally varying complex index of refraction. Often the refractive 
indices are derived from measurements based on in situ dust samples (e.g., Shettle and 
Fenn, 1979; Hess et al., 1998). Optical properties derived from these measurements 
suggest considerably more shortwave absorption by dust than do properties inferred from 
space-based remote sensing techniques (Kaufman et al., 2001; Moulin et al., 2001; Colarco 
et al., 2002; Sinyuk et al., 2003). Retrievals from sun photometer observations also suggest 
very weakly absorbing dust in the shortwave (Kim et al., 2011). Overall, dust absorption as 
represented by the imaginary component of the refractive index may vary by an order of 
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magnitude among these various datasets (Balkanski et al., 2007). Additionally, spectral 
dust optical properties are often derived assuming particle sphericity using Mie theory, 
although it is well known that dust particles are non-spherical. Shape effects are 
understood to be important to retrieval of aerosol properties from spectral, angular, and 
polarized reflectance measurements (e.g., Feng et al., 2009), but the importance to flux 
calculations like those performed in GCM radiative transfer schemes has been thought to be 
small (Mishchenko et al., 1995). Recently, however, improved databases of optical 
properties for non-spherical dust-like aerosols suggest uncertainties in top of atmosphere 
(TOA) fluxes of as much as 30% due to dust shape effects (Yi et al., 2011). While that figure 
is likely on the high end for the effects of non-spherical particles, we are not aware of global 
aerosol modeling studies that have incorporated the effects of dust shape in their radiative 
forcing studies. This provides a motivation for the present study. 

In this study we investigate the impact of dust optical properties on the simulation of the 
dust aerosol lifecycle in the NASA Goddard Earth Observing System (GEOS) atmospheric 
general circulation model. An aerosol and cloud microphysics model has been 
incorporated into the GEOS environment, and parameterizations for dust and sea salt 
aerosol sources and sinks have been included. The tracers in the aerosol module are 
transported online in the GEOS AGCM and are radiatively coupled through the radiative 
transfer scheme, allowing for the inclusion of the aerosol direct radiative effect and the 
associated semi-direct effect on clouds in the model climate. The aerosol indirect effect is 
not considered in this study. A number of simulations are conducted in which we vary the 
dust optical properties both in terms of the particle shape and the spectrally varying index 
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of refraction. In Section 2 we describe the model and the dust optical properties used. In 
Section 3 we evaluate our simulations with observational datasets. In Section 4 we discuss 
the impact of the dust optical properties on the climate of the model. We present 
conclusions in Section 5. 

2. Model Description 

2.1. NASA Goddard Earth Observing System model, version 5 (GEOS-5) 

The GEOS-5 atmospheric general circulation model (AGCM) is one of the main components 
of the GEOS-5 atmospheric data assimilation system and earth system model (Rienecker et 
al., 2008). GEOS-5 is used to provide (i) research quality atmospheric reanalyses for use by 
the scientific community (e.g. the Modern-Era Retrospective Analysis for Research and 
Applications or MERRA, Rienecker et al., 2011), (ii) forward-processed analyses and 
forecasts for use by NASA instrument teams (e.g. Zhu and Gelaro, 2008), and (iii) analyses 
and forecasts for use in various field campaigns (e.g. Bian et al., 2013). It has also been 
used as a tool to study aerosol impacts on weather (e.g. Reale et al., 2011) and climate (e.g. 
Ott et al., 2010; Randles et al., 2013) by inclusion of an online aerosol module (Colarco et 
al., 2010). 

The GEOS-5 AGCM simulates traditional meteorological parameters (wind, temperature, 
pressure, etc.), combining the finite volume dynamical core of Lin (2004) with a column 
physics package. Turbulence is based on Lock et al. (2000). Convection is from the relaxed 
Arakawa-Schubert (RAS) scheme (Moorthi and Suarez, 1992) and incorporates a 
prognostic cloud scheme. Interactions with the land surface are through the catchment 
model of Roster et al. (2000). Radiative transfer is computed using the parameterization of 
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Chou and Suarez (1994; 1999) and Chou et al. (2001), incorporating gaseous absorption 
and scattering and absorption by aerosols in eight shortwave and ten longwave bands. The 
model can be run at various horizontal spatial resolutions. There are 72 vertical layers 
distributed in a hybrid sigma coordinate system that is terrain following near the surface 
and transforming to pressure coordinates near 180 hPa with the model top at 1 Pa 
(~85km). An integrated ocean circulation model was not yet available when this study was 
performed, so sea surface temperatures were prescribed (see below). 


2.2. Community Aerosol and Radiation Model for Atmospheres (CARMA) 

We have coupled the Community Aerosol and Radiation Model for Atmospheres (CARMA, 
Toon et al., 1988; Jensen et al., 1994; Ackerman et al., 1995) to GEOS-5. CARMA is a 
sectional aerosol and cloud microphysics model that has been applied to simulate dust 
(Colarco et al., 2002; 2003; Su et al., 2009), sea salt (Madry et al., 2011; Fan et al., 2011), 
sulfate (English et al., 2011; 2012), and carbonaceous aerosols (Colarco et al., 2004; 
Matichuk et al., 2007; 2008). Although the mechanisms in CARMA permit treatment of 
aerosol growth and nucleation, our setup of CARMA in this study is simple. For dust, we 
partition the aerosol mass across eight logarithmically spaced size bins spanning radius 
from 0.1 - 10 pm. This setup is the same as in Colarco et al. (2003). The dust tracers are 
radiatively coupled to the host model through the GEOS-5 native radiative transfer 
interface. The only CARMA-based microphysical process affecting the dust tracers is 
sedimentation; other processes relating to emissions and deposition are described below. 
Sea salt aerosols are incorporated in a similar fashion (Section 2.5). There is no interaction 


of the dust and sea salt. 
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2.3. Dust Emission and Removal Processes 


Dust emission and removal processes essentially follow from Colarco et al. (2003) and 
Colarco et al. (2010). Dust emissions are based on Ginoux et al. (2001), who provide a map 
of potential dust source locations based on the observed correlation of dust emitting 
regions with large-scale topographic depressions. In this study the topographic map is 
updated and provided at a 1° x 1° latitude-by-longitude spatial resolution (P. Ginoux, 
personal communication) and is remapped in a conservative fashion to the 2° x 2.5° 
horizontal resolution used in our simulations. Dust emissions are computed at each model 
time step and grid box as 

Ei_j s r C ■ S j j u it j{u uj U { ) , u i j Uf ( 1 ) 

where Sij is the grid box i,j potential efficiency to act as a dust source (dimensionless with a 
value 0 - 1 in magnitude) and utj is the model wind speed at 10-meter altitude. u t is the 
particle size-dependent threshold wind speed that must be exceeded in order to emit dust 
(Marticorena and Bergametti, 1995). This threshold wind speed is increased for wet soils 
following Ginoux et al. (2001). s r is the relative fraction of dust particles existing in the soil 
at each of our size bins, similar to Colarco et al. (2003). We assume a globally constant 
distribution for silt-sized particles (1 micron < radius < 10 micron), so s r = 0.25 for each of 
the four size bins in that range. We assume that the mass of clay-sized particles (radius < 1 
micron) is 10% of the silt mass, and the relative partitioning of the four sub-micron size 
bins follows from Tegen and Lads (1996, their Table 1). Cis a global tuning constant used 
to set the overall emissions. In practical terms, since there are no global observations of 
dust emissions, C is tuned to give a reasonable dust loading or dust aerosol optical 
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thickness, and is thus model dependent. In all the experiments described in this study we 
set C = 0.35 pg s' 2 nr 5 . 

Once emitted, dust is transported by advection, as well as turbulent and convective mixing. 
Dust is removed by sedimentation (calculated in the CARMA module), or else by turbulent 
dry deposition and large-scale and convective-scale wet removal, using algorithms that 
were developed for other aerosol packages implemented in GEOS-5 (see Colarco et al., 
2010 ). 

2.4. Dust Optical Properties 

The focus of this study is on the impact of the dust direct radiative effect on the GEOS-5 
atmospheric model. We therefore consider several variations of the dust optical properties 
in our simulations in order to investigate this sensitivity. We approach this by varying the 
assumptions of dust refractive index and particle shape in generating the lookup tables of 
dust optical properties used by the model. Shape is considered for optical purposes only; 
we do not investigate the impact of shape on, for example, sedimentation velocities (e.g., 
Colarco et ah, 2003; Ginoux, 2003). We do not explicitly consider mineralogy or how it 
might vary with source regions, particle size, or during transport. Rather, we prescribe a 
size-independent refractive index that varies only with wavelength. The dust optical 
properties will affect the simulation as far as the radiative feedback on the model, and 
hence the dust lifecycle, but in all other respects the treatment of dust is as described in 
Section 2.3 and is the same for each sensitivity study. 
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The dust refractive indices are drawn from three sources: (1) the commonly used Optical 
Properties for Aerosols and Clouds database (OPAC, Hess et al., 1998), (2) the Shettle and 
Fenn (1979, hereafter referred to as SF) database, and (3) a merger of remote sensing 
based estimates of dust refractive indices in the shortwave with the SF values in the 
longwave (referred to hereafter as OBS). The OBS based refractive indices are from 
Colarco et al. (2002) in the UV (based on simulating the aerosol index from the Total Ozone 
Mapping Spectrometer) and from Kim et al. (2011) in the visible (based on inversions of 
the Aerosol Robotic Network sun photometer observations). A log-power law curve is fit to 
the Colarco et al. (2002) and Kim et al. (2011) refractive indices, joining these smoothly 
with the SF refractive indices at infrared wavelengths. The real and imaginary components 
of the refractive index for each of these data sets are shown in Figure 1. 

The optical properties of spherical dust particles are computed using the Mie theory code 
of Wiscombe (1980). For each of our eight size bins we compute and then integrate the 
optical properties across 10 sub-size bins assuming dM / d(log r) = constant across the size 
bin. The optical properties of non-spherical dust particles are derived from the database of 
Meng et al. (2010), which is a compilation of optical properties for tri-axial ellipsoids based 
on a combination of the solutions derived from Mie theory, the T-Matrix method 
(Mishchenko and Travis, 1998), the discrete dipole approximation (DDA) technique 
(Yurkin and Hoekstra, 2007), and an improved geometric optics method (Yang and Liou, 
1996). The electromagnetic edge effects (Nussenzveig 1992) are included in the extinction 
and absorption cross sections for large size parameters in the geometric optics domain. 
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The nodes in the database are refractive index, size parameter, and two values of particle 
aspect ratio. 

Two shape distributions are considered: a distribution composed primarily of tri-axial 
ellipsoids based on Yi et al. (2011) and a spheroidal distribution from Dubovik et al. 

(2006). Spheroids are a special case of tri-axial ellipsoids and are handled in the Meng et 
al. (2010) database. For a visualization of these shapes see Figure 1 in Meng et al. (2010). 
For each of our size bins we integrate across a sub-size bin distribution similar to what was 
done in the Mie theory calculations, with a second integration across the shape 
distribution. 

Optical properties computed include the mass scattering and extinction efficiencies, 
asymmetry parameter, and the Legendre polynomial moments of the polarized phase 
function. Properties are computed at 61 wavelengths spanning 0.25 - 40 microns (the 
wavelengths considered in the OPAC database). Properties are integrated to the eight 
shortwave and ten longwave bands used in the GEOS-5 radiation codes using a simple 
linear averaging. Ultimately we provide GEOS-5 a lookup table of the relevant optical 
properties (mass extinction efficiency, single scattering albedo, and asymmetry parameter) 
for each size bin and wavelength band. Table 1 summarizes the optical properties used in 
each of our simulations. Figure 2 presents the dust optical properties for each of our 
experiments. Also shown is the deviation of aerosol optical thickness (AOT), single 
scattering albedo (SSA), and asymmetry parameter (g) from a baseline of the OPAC- 


Spheres model assuming the summertime climatological dust aerosol particle size 
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distribution from our OPAC-Spheres model run at Cape Verde (Figure 8, under the main 


Saharan dust plume). The model evolved aerosol particle size distribution has an effective 
radius of 2.24 pm, and the total mass loading has been adjusted to give an AOT of 1 at 550 
nm wavelength for the OPAC-Spheres set of optical properties. Calculations for other dust 
optical tables use the same mass loading and particle size distribution, and so Figure 2 
emphasizes changes in the dust AOT, SSA, and asymmetry parameter due to changes in the 
optical tables. Figure 2 shows results for all 61 wavelengths present in our lookup tables. 
The dust optical properties shown in Figure 2 are summarized for a wavelength of 550 nm 
in Table 1. 


Table 1. Summary of experiments and characteristics. 


Experiment Name 

Refractive Index 

Shape 

AOTssc 

SSA550 

g550 

No Forcing 3 

OPAC b 

Spheres 

1.00 

0.86 

0.75 

OPAC-Spheres 

OPAC 

Spheres 

1.00 

0.86 

0.75 

OPAC-Spheroids 

OPAC 

Spheroids* 5 

1.26 

0.88 

0.77 

OPAC-Ellipsoids 

OPAC 

Ellipsoids f 

1.27 

0.88 

0.77 

SF-Spheres 

Shettle and Fenn c 

Spheres 

1.00 

0.82 

0.76 

SF-Ellipsoids 

Shettle and Fenn 

Ellipsoids 

1.24 

0.84 

0.78 

OBS-Spheres 

Colarco/Kim & Shettle and Fenn d 

Spheres 

0.98 

0.93 

0.74 

OBS-Spheroids 

Colarco/Kim & Shettle and Fenn 

Spheroids 

1.22 

0.92 

0.75 


a In the "No Forcing” run the aerosols are passive tracers. Here the optical properties refer 

to those which are used to compute the diagnostic aerosol optical thickness. 

b OPAC refers to dust refractive indices from Hess et al. (1998). 

c Shettle and Fenn refers to dust refractive indices from Shettle and Fenn (1979). 

d Colarco/Kim & Shettle and Fenn refers to dust refractive indices merged from Colarco et 

al. (2002) and Kim et al. (2011) in the UV/Visible and fit by a log-power law curve to merge 

with the Shettle and Fenn (1979) values in the infrared. 

e The spheroidal shape distribution is from Dubovik et al. (2006). 

fThe ellipsoidal shape distribution is from Yi et al. (2011). 


2.5. Sea salt 

Sea salt aerosol was also simulated simultaneously, but externally mixed with the dust. 
Eight size bins were used here as well, ranging from 0.03 - 10 microns dry radius. The 
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treatment of sea salt is similar to that described in Colarco et al. (2010). The source 
formulation of Gong (2003) was used to determine emissions and the particle size 
distribution. Emissions occur only over open ocean and are a function of the 10-m wind 
speed. Sea salt aerosol is hydrated according to the parameterization of Gerber (1985) as a 
function of size and relative humidity. The hydration of the aerosol affects its fall velocity 
and optical properties. The sea salt aerosols are radiatively coupled to GEOS-5, with 
optical lookup tables generated assuming Mie theory and similar for dust. Refractive 
indices for the dry salt particles are from the OPAC database. The hydrated optical 
properties are computed as a function of relative humidity, with water refractive indices 
from the HITRAN database (Rothman et al., 2009) and assuming a volume average 
refractive index of the salt and water droplet. The sea salt results are nearly identical in all 
of our dust sensitivity simulations, and so we do not discuss them further. 

2.6. Experimental Design 

We have nine potential sensitivity experiments to consider: three sets of refractive indices 
and three shape distributions (spheres, spheroids, and ellipsoids) for each. We have 
completed seven of these simulations that include the radiative effects of dust and sea salt 
aerosols (see Table 1), as well as an additional baseline simulation in which the aerosols 
were not radiatively active in the model (our "No Forcing" run). For each simulation the 
model is initialized with zero aerosol mass concentration beginning in model year 2003. 
The model is run for the period 2003 - 2050, discarding the results from 2003 - 2010 as 
simulation spin-up (i.e., we analyze the results only for the 40 year period 2011 - 2050). 
For all of our simulations we maintain present-day concentrations of greenhouse gases and 
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solar forcing. We use an annually repeating sea surface temperature database based on a 
climatology of sea surface temperatures from El Nino-neutral years (Hurwitz et al., 2011). 
The model is run at 2° x 2.5° horizontal resolution with 72 vertical sigma layers from the 
surface to about 85 km altitude. A total of 16 aerosol tracers are carried (eight for dust and 
eight for sea salt). The dynamical and chemical timesteps were 1800 seconds. 

3. Results 

3.1. Dust Emissions 

Figure 3 shows the climatological global annual mean dust emission magnitude of each of 
our simulations. The mean dust emission magnitude across our simulations is about 2100 
- 2200 Tg yr 1 . There is considerable interannual variability, and the overall range is 
between about 1900 - 2400 Tg yr 1 . Emissions are generally the lowest for the No Forcing 
simulation. Highest emissions are for the simulations using the OPAC refractive indices, 
following by the SF simulations, and finally the OBS simulations, which have similar 
emission magnitude to the No Forcing case. Shape has a minor impact on the simulation 
mean emissions compared to refractive index. The climatological annual mean and 
standard deviation of emissions for each simulation are presented in Table 2. 

Figure 4 shows the climatological, seasonal average dust emissions from our No Forcing 
simulation for December-January-February (DJF), March-April-May (MAM), June-July- 
August (JJA), and September-October-November (SON) over the period 2011 - 2050. 
Strong dust sources are evident across northern Africa, including the sources associated 
with the Bodele Depression and Lake Chad (most active in DJF), the eastern Libyan desert, 
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northern Sudan, and across Mali, Mauritania, Algeria, and Western Sahara (see Prospero et 
al., 2002 for an environmental description of major dust source regions). Other major 
sources appear in northwestern India (most active in JJA), across Saudi Arabia and Iraq, 
and near Lake Eyre in Australia. Less intense dust emission regions include most of 
Australia, northern Mexico and the western United States, across Argentina and Chile, 
southern Africa, and more broadly across Asia and China. The Asian dust emissions appear 
particularly anemic in this model relative to other models (e.g., Ginoux et al., 2004). 
Preliminary studies with higher spatial resolution versions of the GEOS-5 modeling system 
(not shown here) suggest that spatial resolution is a key ingredient in determining realistic 
emissions from Asian dust sources because of the complex interplay of topography and 
surface winds in the Takla Makan and Gobi deserts in China. The coarse resolution of the 
topographic features in our simulations here do not develop the strong downslope wind 
features in these regions that drive real-world dust emissions. Some groups have 
attempted to correct for these deficiencies in coarse spatial resolution climate models by 
imposing distribution functions on the grid box mean wind speeds (e.g., Cakmur et al., 

2004; Su and Toon, 2009). We do not implement such a scheme here. Despite this 
shortcoming over the Asian source regions, we find a higher degree of fidelity between the 
coarse and fine spatial resolution versions of our model for Saharan dust sources, which is 
in any case the globally dominant dust source (Table 2). Accordingly we focus our 
evaluation and analysis on Saharan dust. 

Table 2. Mean aerosol budget diagnostics for our simulations, including emissions, 
atmospheric burden, removal rates, lifetime, and AOT. All quantities are computed globally 
except for North African emissions. 

North 
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Experiment 

Name 

Emissions 
[Tg yr 1 ] 3 

African 
Emissions 
[Tg yr 1 ] 

Burden [Tg] 

kwet [d 1 ] 

kdry [d' 1 ] 

T[d] 

AOT 

No Forcing 

2107 ±93 

1446 ±83 

20.82 ± 0.75 

0.081 

0.194 

3.63 

0.019 

OPAC-Spheres 

2221 ±86 

1556 ± 78 

23.10 ± 0.88 

0.078 

0.187 

3.77 

0.022 

OPAC-Spheroids 

2200 ±83 

1548 ± 69 

23.08 ± 0.88 

0.077 

0.183 

3.85 

0.028 

OPAC-Ellipsoids 

2212 ±80 

1539 ±79 

23.08 ± 0.84 

0.080 

0.190 

3.70 

0.028 

SF-Spheres 

2161 ± 79 

1504 ± 76 

22.66 ± 0.82 

0.078 

0.185 

3.80 

0.021 

SF-Ellipsoids 

2138± 70 

1481 ± 63 

22.50 ± 0.70 

0.076 

0.180 

3.92 

0.028 

OBS-Spheres 

2095 ±76 

1448 ± 69 

20.91 ± 0.65 

0.080 

0.190 

3.70 

0.019 

OBS-Spheroids 

2099 ±87 

1453 ± 78 

20.93 ± 0.79 

0.080 

0.195 

3.65 

0.024 


a The ± indicates the standard deviation. 


3.2. Dust Mass Budget 

In addition to summarizing emissions, Table 2 presents a summary of characteristics of the 
dust mass budget in each of our simulations. The multi-year annual mean atmospheric 
burden of dust follows a similar pattern to the emissions, with the No Forcing and OBS 
simulations having the lowest burden, while the OPAC and SF simulations are the highest. 
There is only about a 5% range in the mean annual emissions between the simulations, 
whereas for the burden the range is closer to about 10% (20.82 Tg for No Forcing versus 
23.10 for OPAC-Spheres). The enhanced loading corresponds to a generally greater 
residence time in the OPAC and SF simulations compared to the No Forcing and OBS 
simulations (e.g., 3.63 days for No Forcing versus 3.92 days for the SF-Ellipsoids 
simulation). Here, residence time is defined as in Colarco et al. (2010) as the burden 
divided by the loss rate. For comparison, the AeroCom-Median model dust lifetime is 4.6 
days (Huneeus et al., 2011) and the GEOS-4 based GOCART dust lifetime was 5.85 days 
(Colarco et al., 2010). Also shown are the dust wet and dry loss frequencies defined as in 
Colarco et al. (2010), which provide an indication of relative importance for loss processes. 
The wet loss rate in our simulations is in the range 0.076 - 0.081 days' 1 , compared to 0.055 


days' 1 in our GEOS-4 simulations (Colarco et al., 2010) and 0.084 days' 1 for the median of 
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the AeroCom models reported in Textor et al. (2006). The dry loss rate is in the range 
0.180 - 0.195 days' 1 here, compared to 0.116 days' 1 for our GEOS-4 simulations and 0.245 
days' 1 for the AeroCom mean. The discussion of the loss rates and lifetimes presented here 
is to put our model in the context of other aerosol models, but there is clearly diversity 
among those various models. Inter-model differences arise from different 
parameterizations of loss processes, as well as model physics, spatial resolution, and 
meteorology. Even with respect to comparison with our GEOS-4 simulations it is important 
to note that the GEOS-4 results came from an older version of our modeling system that 
was driven by meteorological reanalyses and incorporated different model physics, 
whereas here the simulations are from a free-running climate model. 

Figure 5 and Table 3 present a comparison of our climatological annual simulated dust 
surface mass concentration to the multi-year annual mean of measurements at 22 sites. 

The data presented are the same as shown in Huneeus et al. (2011), deriving from the long- 
term measurements managed by the Rosenstiel School of Marine and Atmospheric Sciences 
at the University of Miami (Prospero et al., 1989; Prospero, 1996; Arimoto et al., 1995). 
These sites are located mainly far downwind of dust source regions, and the concentration 
of dust is obtained from high volume air samplers by measuring aluminum concentrations 
after collected filters are ashed in an oven (assuming aluminum content is 8% of total dust 
mass). Figure 5 shows the comparison for our OPAC-Spheres model run, also summarized 
in Table 3 (results for all simulations are in Table SI). Our model is well correlated in the 
multi-year annual mean with the observations, although it generally underestimates the 
mass in the central Pacific sites. Performance appears to be best for sites impacted by 
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Saharan dust, with discrepancies between the model and observations mostly within a 
factor 2-3. For the comparisons to data presented here and elsewhere in this study it is 
important to point out that we are comparing the model climatological fields (whether 
monthly, seasonal, or annual) to long-term averages of the various data sets used, and that 
we are not simulating particular events captured in those data sets. 


Table 3. Comparison of climatological mean dust surface mass concentration from the 
OPAC-Spheres model run to multi-year mean of measurements from Propero et al. (1989), 
Prospero (1996), and Arimoto et al. (1995). 

Site # Site Name Years Lat Lon Elv [m] Observed Model 

fltR m~ 3 ] [HR nr 3 1 

1 Cape Point, South Africa 1992 - 1996 34° 20' S 18°28’E 249 1.78 1.82 

2 Cape Grim, Tasmania 1983 - 1996 40° 40' S 144° 40’ E 94 1.48 2.00 

3 Mawson, Antarctica 1987 - 1996 67° 35' S 62° 30’ E 20 0.10 0.11 

4 Palmer Station, Antarctica 1990- 1996 64° 46' S 64°02’W 10 0.03 0.13 

5 Yate, New Caledonia 1983 - 1985 22° 08' S 167° 00' E 0 0.17 0.22 

6 Funafuti, Tuvalu 1983 - 1987 08° 30' S 179° 11’ W 0 0.20 0.01 

7 Nauru 1983 - 1987 00° 31' S 166° 56' E 7 0.10 0.01 

8 Norfolk Island 1983 - 1997 29°04’S 167° 58' E 0 0.84 0.39 

9 Rarotonga, Cook Islands 1983 - 1994 21° 30' S 159° 45’ W 0 0.10 0.11 

10 American Samoa 1983 - 1996 14° 15' S 170° 34’ W 42 0.16 0.03 

11 Midway Island 1981 - 1997 28° 13' N 177° 20' W 0 0.72 1.20 

12 Cheju, South Korea 1991 - 1995 33°31’N 126° 28' E 70 14.14 2.91 

13 Hedo, Japan 1991,1994 26° 55' N 128° 15' E 0 8.37 1.54 

14 Fanning Island 1981 - 1986 03° 53’ N 159° 19’ W 4 0.10 0.02 

15 Enewetak Atoll 1981 - 1987 11° 19' N 162° 19' E 0 0.24 0.12 

16 Ragged Point, Barbados 1984- 1998 13° 10’ N 59°25’W 41 16.15 27.18 

17 Izana, Tenerife 1987 - 1998 28° 18' N 16°30’W 2367 30.18 74.45 

18 Bermuda (East & West) 1989 - 1998 32° 16’ N 64°52’W 51 3.36 6.93 

19 Mace Head, Ireland 1988 -1994 53° 19' N 09°50’W 5 1.01 0.91 

20 RSMAS - University of Miami 1989 - 1998 25°45’N 80°15’W 5 4.59 11.38 

21 Rukomechi, Zimbabwe Not Known 16° 00’ S 29°30’E Unknown 10.53 0.21 

22 Jabirun, N. Australia Not Known 12°41’S 139° 53’ E Unknown 4.03 2.98 


Comparisons of our model to multi-year estimates of dust deposition from the OPAC- 
Spheres model run are shown in Figure 6, with site locations and details for all of our 
model runs presented in Tables S2 and S3. Again, the data are as presented and screened 
in Huneeus et al. (2011) and are derived from measurements taken during the Sea/Air 
Exchange field campaign (SEAREX, Prospero et al., 1989, with data as presented in Ginoux 
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et al., 2001), measurements presented in Mahowald et al. (2009), and the Dust Indicators 


and Records in Terrestrial and Marine Paleoenvironments dataset (DIRTMAP, Tegen et al., 
2002; Kohfeld and Harrison, 2001). The model is generally underestimating dust 
deposition, particularly in the remote Pacific locations, and is overestimating in the remote 
regions near Antarctica, but is otherwise generally within a factor of 10 of the observations. 


3.3. DustAOT 

In Figure 7 we show the climatological, seasonal average dust AOT (at 550 nm) from our 
No Forcing simulation. Consistent with Figure 4 we see high dust AOT associated with the 
Bodele Depression in northern Africa peaking in DJF, the strong dust sources in western 
India peaking in JJA, and the peak of the Asian dust emissions during MAM. The most 
evident feature remote from the source regions is the Saharan dust plume over the North 
Atlantic Ocean. This plume moves north and south with the seasonal progression of the 
inter-tropical convergence zone (ITCZ), with the peak of both its northward and westward 
transport occurring during JJA. 

Figure 8 shows a comparison of our model 550 nm dust AOT over the northern Atlantic 
Ocean to observations retrieved by the Multi-angle Imaging Spectroradiometer (MISR) 
flying onboard the NASA Terra spacecraft. MISR measures reflected radiation in four 
visible channels from nine different cameras viewing the same scene at different angles as 
the spacecraft flies over. These 36 pieces of information permit cloud clearing and 
selection of an appropriate lookup table of aerosol radiative properties to retrieve AOT, 
among other aerosol characterization capabilities (Kahn et al., 2005). Since our model is 
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not simulating the actual period of MISR observations we compare here our JJA model 
climatological dust AOT to the JJA climatology of MISR AOT for the period 2000 - 2011. For 
contrast with Figure 7, we show in Figure 8 the simulated dust AOT from both our No 
Forcing and our OPAC-Spheres model runs. Comparing these we clearly see the higher 
AOT in the OPAC-Spheres run consistent with higher dust emissions, as well as a further 
westward transport of the dust plume over the Caribbean. In comparison to MISR, the 
OPAC-Spheres model run has a similar north-south placement of the dust plume, although 
the model AOT appears somewhat higher than the observations over the source regions. 
The MISR observations are additionally sensitive to other aerosols besides dust, as is 
evident in the pollution plume along the U.S. east coast and the features associated with 
biomass burning in South America and southern Africa. 

Figure 9 shows the climatological JJA simulated 550 nm dust AOT for each of our 
simulations at four different longitudes west of northern Africa. Also shown is the 2002 - 
2011 JJA climatology of the 550 nm AOT from the Moderate Resolution Imaging 
Spectroradiometer (MODIS) onboard the NASA Aqua spacecraft. MODIS retrieves over 
ocean AOT in cloud-free conditions from measurements of reflected radiation at six visible 
channels (Tanre et al., 1997). Similar to the MISR observations shown, we are not 
attempting to separate the dust AOT from the total AOT in the MODIS observations, but 
rather focusing on what the observations tell us qualitatively about our simulations. Near 
the source region (the slice at 20° W) the different simulations all have a peak in the dust 
AOT at about 20° N, consistent with the placement of the peak in the MODIS AOT. AOT is 
highest for the OPAC based simulations using non-spherical optics, consistent with the 
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higher emissions in these simulations (Figure 3) and the enhanced mass extinction 
efficiency of the non-spherical optics (Figure 2). The north-south placement of the dust 
plume is similar for all of the model simulations, and consistent with the peak in the MODIS 
AOT out to 60° W. By 80° W this similarity with the MODIS observations has broken down. 
As the plume moves to the west the No Forcing and OBS based simulations (i.e., the least 
absorbing simulations) peak in AOT more toward the southern edge of the dust plume, 
while the OPAC and SF based simulations maintain a more uniform north-south 
distribution and higher AOT magnitude. 

Figure 10 shows the climatological JJA simulated 550 nm dust AOT for each of our 
simulations across the North Atlantic Ocean averaged in the latitude band 10° - 30° N. For 
comparison we also show the JJA MODIS Aqua over ocean total AOT averaged over the 
period 2002 - 2011. In Figure 10a the results are plotted exactly according to the model 
output, and the difference in the model simulations due to optical assumptions and 
differences in emissions is clearly evident at the eastern longitudes, closest to the source 
regions. In Figure 10b the model and MODIS AOT are normalized to the same value of AOT 
at the longitude of Cape Verde (off the west coast of northern Africa, near the source 
region, see Figure 8) and the results are plotted on a log-y axis. This emphasizes the 
differences in the long-range transport of the dust from northern Africa, across the Atlantic, 
and into the Caribbean, showing that in a relative sense dust is transported further west in 
the OPAC and SF based simulations, while the OBS and No Forcing simulations are more 
similar with the least westward export of the dust plume. Also evident is that the relative 
westward transport of the dust plume is largely insensitive to the shape assumptions, with 


473 

474 

475 

476 

477 

478 

479 

480 

481 

482 

483 

484 

485 

486 

487 

488 

489 

490 

491 

492 

493 

494 

495 


the model results most clearly clustered according to refractive index in Figure 10b. The 
MODIS observations lend some context, showing the diminishment of AOT as the plume 
propagates west of Africa, although clearly the model results fall off more rapidly than the 
satellite observations suggest. This could be due to too aggressive removal of aerosol in the 
model. Alternatively, the MODIS observations of course also include aerosols other than 
dust, which become an important component of the total AOT in the westernmost 
longitudes (see Figure 8b). 

Figure 11 shows a comparison of the climatological monthly mean 550 nm dust AOT from 
each of our simulations compared to climatological monthly averaged observations from 
the Aerosol Robotic Network (AERONET) of sun/sky scanning photometers (Holben et al., 
1998). AERONET observes attenuation of the direct solar beam in measuring AOT at a 
number of wavelengths, here interpolated to 550 nm wavelength. The AERONET site 
locations chosen are shown in Figure 8 and described in Table 4, which also shows the 
years and total number of monthly mean AOT values being used in Figure 11. Near the 
source region (Santa Cruz, Dahkla, Cape Verde, and Tamanrasset) the model has a similar 
seasonal cycle to the observed AOT, although the model tends to have a higher magnitude 
AOT relative to the observations, especially at Dahkla and Santa Cruz. At Tamanrasset the 
simulated AOT is more comparable to the magnitude in the observations, although the 
seasonal cycle in the model is distinctly bimodal in all simulations while it is relatively flat 
between March and September in the observations. On the other hand, the seasonal cycle 
and magnitude are best simulated at Cape Verde, which is directly below the main 


summertime transport pathway. 
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Table 4. Description of AERONET sites used in this study. Included are the responsible 
principle investigator (PI) at each site, the years active, the latitude (Lat), longitude (Lon), 
and elevation (Elv). Also indicated are the number of monthly mean AOT values used in 
constructing Figure 11 (n ao t) and the number of JJA monthly mean inversions used in 


constructing the seasonal particle size distributions shown in Figure 12 (nmv). 


Site 

PI 

Years 

Lat 

Lon 

Elv [m] 

llaot 

llinv 

Ragged Point 

Joe Prospero 

2007 -2011 

13° 09' N 

59° 25' W 

40 

52 

n 

La Parguera 

Brent Holben 

2000 -2011 

17° 58' N 

67° 02' W 

12 

110 

24 

Bermuda 

Brent Holben 

1996-2002 

32° 22 N 

64° 41 W 

10 

64 

10 

Cape Verde 

Didier Tanre 

1994-2011 

16° 43' N 

22° 56' W 

60 

179 

39 

Santa Cruz 

Emilio Cueva-Agullo 

2005-2011 

28° 28' N 

16° 14 W 

52 

72 

17 

Dahkla 

Hammad Benchekroun 

2002 -2003 

23° 43' N 

15° 56' W 

12 

22 

6 

Tamanrasset 

Emilio Cueva-Agullo 

2006-2009 

22° 47' N 

5° 31' E 

1377 

24 

6 


Figure 11 also shows several sites impacted on the western edge of the dust plume 
(Bermuda, La Parguera, and Ragged Point). Bermuda is only marginally impacted by dust, 
as evidenced by the relatively high Angstrom parameters at that site (not shown). There is 
a clear difference here again between the No Forcing and OBS based simulations on the one 
hand and the OPAC and SF simulations on the other, showing the more pronounced 
northward transport in the OPAC and SF simulations while the No Forcing and OBS based 
simulations have a pronounced minimum in AOT during the summer months. At La 
Parguera the summertime Angstrom parameter is generally low, indicative of dust, and the 
seasonal peak in dust AOT magnitude in June and July is well representing in the model, 
although the model appears to shift the peak toward the earlier part of the season. Again, 
the distinction between the more absorbing SF and OPAC simulations versus the No 
Forcing and OBS simulations is clear, with much lower AOT in the latter set. This remains 
the case at Ragged Point, the furthest south site considered, where again the seasonal cycle 
and peak AOT magnitude are well represented by the model. At all three sites we compare 


the ratio of the June climatological simulated dust AOT between the SF-Spheroids and No 
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Forcing simulations. The ratio is highest Bermuda (1.78), intermediate at La Parguera 
(1.59), lowest at Ragged Point (1.34). This indicates a relative preference for southward 
transport of the dust plume in the No Forcing model, which is also the case for the OBS 
based models. Incidentally, this ratio is lower than the Ragged Point value at all the near- 
source sites except Tamanrasset (where it is 1.44). 

3.4. Dust Particle Size Distribution 

Figure 12 shows a comparison of the volume dust particle size distribution simulated in the 
model to the retrieved particle size distribution from the AERONET measurements. Under 
high AOT (> 0.4) and homogeneous sky conditions AERONET can invert column integrated 
particle size distribution from its almucantar scan (Dubovik et al., 2000). For each of the 
sites shown we aggregated monthly mean particle size distributions from the AERONET 
dataset during JJA (the number of individual months used at each site is indicated in Table 
4). For comparison we show the climatological JJA column integrated particle size 
distribution from the model. Note that the AERONET retrievals all contain a fine mode 
component to the total particle size distribution (radius < 0.5 um) that is absent in the 
model. This fine mode presumably corresponds to aerosols not simulated in the model. 
Near the source regions (Cape Verde and Tamanrasset) there is very little difference in the 
particle size distribution among the model runs. The model tends to peak in volume at 
radius ~ 1.5 um while the AERONET retrievals put the volume peak between about 2-3 
um radius. Farther from the source regions (La Parguera and Ragged Point) the peak in the 
simulated particle size distribution volume occurs in the same size bin as near the sources. 
However, larger particles, preferentially removed by sedimentation in the model, are 
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relatively depleted compared to the retrievals at these distant sites. There are as well 
larger differences between the model simulations at these sites, owing to greater difference 
in the total loading at these sites. 

3.5. Dust Vertical Profiles 

Figure 13 shows a composite of the climatological JJA dust vertical profile from three of our 
model simulations at four longitudes west of northern Africa. The dust vertical profile is 
presented as the profile of the total dust mass mixing ratio. Also shown is the JJA 2011 
aerosol attenuated backscatter profile measured by the Cloud-Aerosol Lidar with 
Orthogonal Polarization (CALIOP) flying aboard the NASA CALIPSO spacecraft. CALIPSO 
flies in a polar orbit in formation with the NASA Aqua spacecraft as part of the so-called "A- 
Train" of satellites. CALIOP profiles the atmosphere with a two-channel lidar system, 
measuring attenuated backscatter at 532 and 1064 nm and depolarization at 532 nm 
(Winker et al., 2010). The comparison of the model vertical profile to the CALIOP 
observations should be understood to be qualitative, as only a single season of the CALIOP 
observations are presented, and of course the CALIOP backscatter observations are 
sensitive to the presence of clouds and other aerosols besides dust. Still, the dust features 
seen in the CALIOP observations are fairly typical of the observed profiles, showing an 
elevated aerosol layer near the African coast (20° W) between about 850 - 550 hPa altitude 
and between about 10° - 25° N latitude. This is similar to the vertical placement of the 
aerosol plume in the model simulations shown, which peaks at about 20° N latitude and 
between about 850 - 700 hPa altitude. As the aerosol moves west the elevated dust plume 
is still evident in the CALIOP observations at 40° W at nearly the same latitudes, although 
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in the model the plume has moved south to a centroid of around 15° N latitude. Further 
west the aerosol layer is less distinct in the CALIOP observations and essentially merges 
with boundary layer aerosols. The simulated profile descends further toward the west, and 
it becomes evident that the dust retains its elevation more clearly in the OPAC-Spheroids 
simulation while in the OBS-Spheroids and No Forcing simulations it is at a lower altitude. 

Figure 14 shows the northern Atlantic JJA climatological dust median altitude (the altitude 
at which half the dust column mass is above and below) for each of our simulations. This 
approach normalizes differences in the absolute column loading between the simulations. 
Over the dust source region, where aerosol mass is well mixed by the deep boundary layer, 
there is little difference between the simulations. When the dust is advected over the 
ocean, however, the differences in the simulations are immediately apparent, with the 
OPAC and SF simulations maintaining a much higher dust median altitude than either the 
OBS or No Forcing simulations. The difference in the median altitude is about 100 hPa 
between the SF and No Forcing simulations at about 60° W. These results are consistent 
with what is shown in Figure 13 and suggest an explanation for the higher AOT maintained 
by the SF and OPAC simulations across the Atlantic as shown in Figure 10: the relatively 
higher in altitude the dust mass remains during transport, the longer its atmospheric 
residence time and hence the further west it may be transported. 

4. Discussion 


4.1. Summary of model results 
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In Section 3 we presented our simulations of the dust aerosol lifecycle using the CARMA 
aerosol module run online in GEOS-5. Model results were shown for the 40-year period 
2011 - 2050 under conditions of a repeating annual cycle of sea surface temperatures 
taken from an El Nino neutral climatology. A baseline simulation had aerosol radiative 
interaction turned off in the model. In all other simulations dust radiative interaction with 
the AGCM was considered, with assumptions about dust optical refractive index and shape 
distribution varied between the simulations. 

For all simulations the climatological global, annual mean dust emissions were in the range 
of 2095 - 2221 Tg yr 1 (Table 2). North African dust emissions were in the range 1446 - 
1556 Tg yr 1 . The difference between the high and low mean emissions are thus about 6%. 
Our emissions are similar in magnitude to the multi-year simulations presented in the 
GEOS-4 based version of our modeling system that employed essentially the same dust 
mobilization scheme (Colarco et al., 2010). There the multi-year average dust emissions 
over the same particle size range were 1970 Tg yr 1 . The global, annual mean dust aerosol 
burden here is in the range 20.8 - 23.1 Tg. This is less than the burden of 31.6 Tg found in 
Colarco et al. (2010), and the aerosol residence time is correspondingly found to be shorter 
in the present study (3.63 - 3.92 days) than in the previous one (5.85 days). There are 
important differences between the GEOS-5 based model and the GEOS-4 system: the 
physical parameterizations for turbulence, convection, and precipitation are entirely 
different, which have significant impacts on the aerosol transport and removal, and the 
GEOS-4 simulations were run with assimilated meteorology, whereas here we are running 
climate simulations. The shorter dust aerosol atmospheric residence time and lower dust 
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burden here are more consistent with the mean of models participating in the Aerosol 
Comparisons between Models and Observations (AeroCom) project (19.2 Tg and 4.22 days, 
respectively, see Textor et al., 2006; see also the updated values of 15.8 Tg and 4.6 days, 
respectively, in Huneeus et ah, 2011). Climatological annual mean dust surface mass 
concentrations were within a factor of 3 of multi-year observations for sites impacted 
primarily by long-range transport of Saharan dust. At remote sites in the Pacific the model 
tended toward underestimating the observed values (Figure 5). The model tends to 
underestimate deposition at sites impacted by Saharan dust (Figure 6), albeit mostly within 
a factor of 3 - 5, but overestimates deposition at Antarctic sites, although the fluxes are 
quite low. 

Despite similarities in the dust emissions, burden, and lifetime, there is a wide range in the 
value of the global, annual mean dust AOT in our simulations (range is 0.019 - 0.028). The 
higher dust AOT values are associated with the model runs that incorporated non-spherical 
dust optics, owing to the higher visible wavelength mass extinction efficiency (i.e., AOT per 
unit mass) in the non-spherical cases (see Figure 2). The various simulations show similar 
AOT spatial and temporal variability as satellite and ground-based observations (see 
Figures 8-11), although there a tendency of the model AOT to be somewhat higher in 
magnitude than the observations, particularly near the Saharan dust sources. Differences 
in the particle size distribution are manifest mainly as differences in the loading (Figure 
12) with a generally similar modal radius and width of the size distribution among the 
various simulations. The vertical distribution of dust downwind of the Sahara differs 
significantly among the simulations, with similar altitudes of the dust plume near the 
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continental sources (Figure 13) but a clear differentiation during westward transport 
(Figure 14). Median dust altitude is lowest for No Forcing simulation and only slightly 
higher for the OBS simulations, while the SF and OPAC simulations instead have distinctly 
higher dust median altitude. 

Consistent with the vertical distribution discussed above, the east-west gradient in dust 
AOT across the Atlantic is more moderate in the OPAC and SF simulations than for the OBS 
and No Forcing runs (Figure 10b). This gradient is not as moderate as the MODIS 
observations, suggesting that removal mechanisms in the model may be somewhat 
overestimated. Previous studies with our model have investigated the east-west gradient 
of dust AOT across the Atlantic and determined wet removal was too aggressive 
(Nowottnick et al., 2011). While we adopted Nowottnick et al.’s (2011) suggested changes 
to the scavenging parameters in these simulations it seems this issue is persistent. There 
are differences between this study and the Nowottnick et al. (2011): their results were 
from a single seasonal simulation of a higher resolution run of the model driven by 
assimilated meteorology. Nevertheless, the more moderate east-west gradient in AOT, the 
higher transport altitude, and the more northward shift of the dust plume in the OPAC and 
SF simulations suggest they agree better with the AERONET observations (Figure 11) than 
do the OBS and No Forcing simulations, and thus provide the most realistic simulation. 

4.2. Dust Radiative Forcing 

We extend the analysis of the simulated dust AOT to encompass the dust direct radiative 
effect on the climate in GEOS-5. Our instantaneous forcing follows the presentation in 
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Randles et al. (2013), and is calculated as the net (positive down) flux change at the top of 


atmosphere (TOA) with and without aerosols, holding the atmospheric state fixed between 
two separate calls to the radiative transfer module. Defined this way, a positive TOA 
forcing indicates the addition of energy to the climate system (i.e., a radiative warming 
effect), while a negative forcing indicates a net energy loss (i.e., a radiative cooling effect). 
Table 5 summarizes the climatological JJA global clear-sky aerosol radiative effect in our 
simulations (the global, annual mean climatological values are presented in Table S4). 


Table 5. Global, climatological JJA clear-sky aerosol radiative effect over ocean (land) [W 
nr 2 ]. 



SW TOA 

SW ATM 

SW SFC 

LW TOA 

LW ATM 

LW SFC 

TOT ATM 

OPAC-Spheres 

-1.62 [-0.01) 

1.04 [4.44) 

-2.66 [-4.45) 

0.22 [0.55) 

-1.08 [-1.49) 

1.30 [2.04) 

-0.04 (2.95) 

OPAC-Spheroids 

-1.78 [0.07) 

1.26 [5.29) 

-3.04 [-5.22) 

0.22 [0.53) 

-1.07 [-1.46) 

1.29 [1.98) 

0.20 [3.83) 

OPAC-Ellipsoids 

-1.66 [-0.03) 

1.11 [4.71) 

-2.76 [-4.74) 

0.22 [0.53) 

-1.07 [-1.47) 

1.29 (1.99) 

0.04 [3.24) 

SF-Spheres 

-1.46 [0.62) 

1.20 [5.22) 

-2.66 [-4.60) 

0.17 [0.25) 

-0.96 [-0.80) 

1.14 [1.05) 

0.23 (4.41) 

SF-Ellipsoids 

-1.49 [0.57) 

1.26 [5.42) 

-2.75 [-4.84) 

0.17 [0.22) 

-0.95 [-0.74) 

1.12 (0.96) 

0.31 (4.67) 

OBS-Spheres 

-1.49 [-0.60) 

0.49 [2.12) 

-1.98 [-2.72) 

0.16 [0.22) 

-0.96 [-0.75) 

1.13 [0.97) 

-0.48 [1.37) 

OBS-Spheroids 

-1.50 f-0.54) 

0.58 [2.56) 

-2.09 [-3. 10) 

0.16 [0.20) 

-0.95 f-0.70) 

1.11 [0.90) 

-0.37 (1.86) 


Table 5 indicates a net clear-sky TOA SW cooling over the oceans for all simulations. Over 
land the picture is somewhat more complicated, depending on the dust optical model 
chosen. The weakly absorbing OBS models also exert a SW TOA clear-sky cooling over 
land, while the more strongly absorbing SF models exert a warming. The OPAC models 
have nearly zero net TOA effect over land. Figure 15 shows the climatological, JJA TOA 
"forcing efficiency" (forcing per unit AOT) for several of our model runs over northern 
Africa and the North Atlantic Ocean. Note that we show forcing efficiency in Figure 15 to 
facilitate comparison of the simulations, but that absolute forcing values are shown in 
Tables 5 and S4. So, for example, we show the forcing efficiencies for our SF-Spheres and 


OPAC-Spheres cases, as well as the difference between them (SF-Spheres minus OPAC- 
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Spheres). For both simulations there is the negative forcing efficiency over the ocean and 
positive forcing efficiency over the dust source regions, and the difference plot is consistent 
with the values in Table 5. 

Also shown in Figure 15 is the shape effect, presented as the difference of the OPAC- 
Ellipsoids and OPAC-Spheroids models relative to the OPAC-Spheres baseline. The non- 
spherical optical models are relatively cooling at TOA over the source region compared to 
the spherical model, but are somewhat less cooling over the oceans or land remote from 
the source regions (positive differences in Figure 15). The presentation in Figure 15 is not 
obviously consistent with the values in Table 5, but we remind that we are presenting 
global mean absolute forcing in Table 5, while it is forcing efficiency shown in Figure 15. 

Finally, in Figure 15 we also show the difference of the OBS-Spheres model compared to 
the OPAC-Spheres baseline. Here the OBS-Spheres model is shown to more efficiently 
cooling over both and land and ocean, although there is a warming feature over the 
southern portion of the dust plume, consistent with the relatively more southward 
transport of the dust in the OBS-Spheres model run (see Figure9). 

The dust clear-sky longwave TOA effect is warming for all simulations, as shown in Table 5. 
The global average warming over ocean is relatively small and similar for all simulations, in 
the range of 0.16 - 0.22 W nr 2 . Over land the effect is larger and more varied, between 0.20 
- 0.55 W nr 2 . The strongest longwave effect is for the OPAC based simulations, while the SF 
and OBS simulations are more similar. Recall that the OBS and SF simulations used the 
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same refractive indices in the longwave. Figure 2 shows the relatively higher longwave 
AOT and lower single scattering albedo for the OPAC refractive indices compared to the 
other simulations, hence the stronger longwave effect. 

Surface effects are consistent among the simulations, with an overall cooling in the 
shortwave and warming in the longwave, and a greater magnitude for both over land than 
ocean. 

The atmospheric forcing is the difference between the TOA and surface terms in Table 5, 
and is the amount of energy added or subtracted from the atmosphere due to the dust 
radiative effects. Here the distinction between the different refractive index sets becomes 
clear, with the OBS simulations having total forcing in the range -0.37 to -0.48 W nr 2 over 
ocean and 1.37 - 1.86 W nr 2 over land. These can be compared to the relatively warmer 
atmospheric forcing in the OPAC (-0.04 to 0.20 W m-2 over ocean, 2.95 - 3.83 W nr 2 over 
land) and the SF (0.04 - 0.23 W nr 2 over ocean, 4.41 - 4.67 W nr 2 over land) simulations. 
The SF simulations have overall the greatest atmospheric direct effect, owing to their 
weaker longwave atmospheric cooling. 

Differences in the direct radiative effect among the simulations due to particle shape are 
not clear in the longwave, although there is some consistency among them in the 
shortwave. The spheroidal simulations have relatively stronger shortwave TOA and 
surface cooling over and ocean compared to the associated spherical simulations, are 
slightly more warming over land, and have relatively greater atmospheric absorption. The 
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ellipsoidal simulations also indicate a stronger atmospheric forcing than the spherical 
simulations, albeit a smaller effect than for the spheroidal models. The ellipsoids also have 
slightly higher TOA cooling over ocean and in contrast with the spheroids are also 
relatively cooling over land. 

4.3. Dynamical Impact of Dust Heating 

The radiative forcing of dust feeds back onto the dynamics of the GEOS-5 model. An 
obvious feedback is on the cycle of dust mobilization itself. For example, as Figure 3 and 
Table 2 illustrates there is an impact of the dust optical properties on both the magnitude 
and interannual variability in dust emissions. Emissions were found to be on average 
lowest for the No Forcing simulation and then increasing in magnitude for the OBS, SF, and 
then OPAC based simulations. This result is in contrast to previous studies that found 
decreasing dust emissions for more absorbing dust aerosols (e.g., Perlwitz and Miller, 

2001; Miller et al., 2004). The emissions magnitude was more related to the refractive 
index choice than particle shape, a grouping made clearer in the total atmospheric forcing 
column of Table 5 (TOT ATM = SW ATM + LW ATM). Dust atmospheric forcing is highest 
for the SF simulations and lowest for the OBS simulations, with OPAC slightly less than the 
SF simulations. 

Figure 16 shows the JJA climatological difference in the 10-m winds between the OPAC- 
Spheres and No Forcing model runs, as well as for the OBS-Spheres and No Forcing model 
runs. The general pattern is for an enhancement in the climatological mean wind field near 
the southern edge of the Sahara for the simulations that include dust radiative forcing 
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relative to the No Forcing simulation. The band of higher mean wind speed extends across 
the western half of the continent between about 15° - 20° N, and is partially offset by 
relatively weaker winds in the dust forcing simulations to the north. The wind vectors 
shown in Figure 16 indicate this enhanced wind speed across the southern Sahara is 
generally an onshore flow pattern, resulting from enhanced heating over land due to dust 
relative to the fixed (cooler) sea surface temperature immediately offshore. This is 
consistent with the previous study of Lau et al. (2009). 

Also shown in Figure 16 are the JJA climatological differences in dust emissions between 
the two simulation pairs. Emissions are the convolution of the surface wind speeds with 
the dust source function (the thin contours in Figure 16, see also Equation 1). The region of 
enhanced surface winds are correlated with relatively enhanced dust emissions in the dust 
forcing simulations. This relationship is not perfect as the mean wind speeds do not tell the 
whole story: although the mean wind speed difference is small or even negative to the 
north of this band there are still generally enhanced emissions in this region, especially in 
the OPAC-Spheres simulation. Because of the near cubic dependence of the emissions on 
the surface wind speed, emissions are sensitive to the distribution of wind speeds rather 
than simply to the mean. 

Figure 17 presents a complementary look at the dynamics of this system. Here we show 
the JJA climatological mean difference in the vertical profile of atmospheric temperature for 
the same two pairs of simulations. The profile is the zonal mean between 15° W and 0° 
longitude (the westernmost portion of the dust source region). Also shown is the 
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difference in the zonal mean wind speed between the forcing and No Forcing simulations 
(thin solid and dashed contours) and the isosurface of dust concentrations exceeding 100 
pg nr 3 in dust concentration (thick contour line). For both pairs of simulations the 
presence of dust acts to heat the atmosphere between about 800 - 600 hPa altitude. 
Attenuation of incident shortwave radiation leads to a cooling under the dust plume, while 
longwave cooling at the top of the dust layer leads to a cooling between about 500 - 400 
hPa altitude. This temperature signal is present in both pairs of simulations, although the 
signal has much higher amplitude in the OPAC-Spheres simulation, which has considerably 
greater atmospheric forcing and surface cooling than the OBS-Spheres simulation (Table 
5). The modified temperature profile in the dust forcing simulations alters the zonal winds, 
adding a westerly component at the surface but strengthening the easterly jet aloft (by 
about 2.5 m s _1 for the OPAC simulations), with the peak increase in the easterly winds at 
about 550 hPa altitude on the southern side of the dust plume. This result is consistent 
with the findings in Kim et al. (2010), who used a GCM with climatological aerosols to 
identify a strengthening of the African easterly jet south of Saharan air layer resulting from 
thermal wind balance with the addition of dust heating. 

Figure 18 puts the dust heating effect into its global context. Here we show the JJA 
climatological zonal winds and mean meridional circulation (MMC) for our No Forcing 
simulation (contours), as well as the difference (shading) between our dust forcing 
simulations (again, OPAC-Spheres and OBS-Spheres) and the No Forcing model run. For 
the zonal mean wind, the signals seen in Figure 18 are similar to those shown for natural 


aerosols in Randles et al. (2013), who used a somewhat older version of our aerosol 
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modeling system (their Figure 7h). For both of our dust forcing simulations shown here, 
there is the enhancement in the mid-tropospheric easterly wind speed in the main Saharan 
dust plume (between 700 and 500 hPa and 0° to 30° N, similar to Figure 17). There is as 
well a strengthening and slight poleward shift of the southern hemisphere jet stream (at 
200 hPa and 30° S). Both of these signals are stronger in magnitude in the OPAC-Spheres 
simulation than the OBS-Spheres. There is also an enhancement of the westerlies at about 
60° S, which is similar in magnitude and structure for both pairs of our simulations. As this 
is a mostly dust-free environment this enhancement is likely due to the sea salt aerosols, 
which were very similar in all simulations (but not radiatively coupled in the No Forcing 
simulation). Note that the southern hemisphere enhancement in the westerlies is stronger 
in Randles et al. (2013), who also had a significant high bias in sea salt AOT that has largely 
been corrected in our model. 

Figure 18 also shows the mean meridional circulation, a representation of the meridional- 
vertical mass stream function. The contours in Figure 18 show the sign of the MMC for the 
No Forcing simulation, with the arrows giving the sense and qualitative magnitude of the 
associated vertical motion. The effect of the aerosol forcing is to strengthen both the 
northern and southern branches of the Hadley cell, enhancing both the upward motion on 
the equatorward side of the cells and the downward motion on the poleward side. The 
perturbations suggest a slightly northward shift in the cell's positions, supported also by a 
northward shift of the ITCZ evident in the modeled precipitation fields (not shown), similar 
to Wilcox et al. (2010). Again, the strengthening and the shift are more pronounced for the 
OPAC-Spheres simulation than the OBS-Spheres case. 
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The results for the MMC are somewhat at odds with Randles et al. (2013), who diagnosed a 
summertime weakening of the MMC due to natural aerosols, but there are some significant 
differences between their simulations and ours. For one, they had a much higher baseline 
forcing due to high sea salt AOT (noted above). Additionally, they had annually varying 
SSTs, while we are using a fixed, repeating cycle of SSTs. Finally, their shortwave dust 
forcing was weak compared to the OPAC-Spheres model, more like what is shown in the 
OBS-Spheres, which has a weaker impact overall. On the other hand, the Randles et al. 
(2013) were broadly consistent with Allen and Sherwood (2010) who showed a similar 
weakening of the JJA Hadley cell due to dust and sea salt for climate model simulations with 
imposed (not dynamic) aerosol concentrations and fixed SSTs. Our results for the MMC, by 
contrast, are more like what Randles et al. (2013) and Allen and Sherwood (2010) show for 
anthropogenic aerosols, strengthening the Hadley cells. The anthropogenic aerosols in 
those studies are typically quite absorbing in the shortwave, in contrast to our OBS-Spheres 
simulations but more like our OPAC-Spheres model run. 

5. Conclusions 

The discussion above suggests several general conclusions of our study. First, we find that 
higher Saharan dust emissions are associated with increasing the radiative impact of the 
dust. Table 5 shows greatest overall atmospheric heating for the SF based simulations, 
followed by the OPAC and the OBS based simulations. Dust emissions are similar, highest 
for the OPAC simulations, lower for the SF runs, and lowest for the OBS and No Forcing 
simulations (Figure 3 and Table 2). The higher dust emissions are the result of 
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enhancement in the surface wind speeds over important dust source regions (Figure 16). 
These results are opposite to studies by Perlwitz and Miller (2001) and Miller et al. (2004), 
where the highest magnitude dust emissions were found for simulations in which dust was 
a passive tracer and the lowest overall emissions were found for the most absorbing dust 
optical models. Those studies were carried out in a relatively coarser spatial resolution 
AGCM (4° x 5° latitude by longitude with 12 vertical layers), and the Miller et al. (2004) 
study included a mixed-layer ocean model, so there are important differences in the overall 
experiment construction. 

Second, we find that dust is transported more efficiently and at higher altitudes in our 
model the more absorbing the dust optical properties. Over the northern Atlantic Ocean 
dust is transported further to the west and at a more northern latitude for our OPAC based 
simulations when compared to the No Forcing simulation (Figures 9 and 10). Relatively 
furthest westward transport of the summertime Saharan dust plume is in our SF based 
simulations, while the weakly absorbing OBS based simulations are most like the No 
Forcing case. Dust median altitude is highest for the SF and OPAC simulations, and lowest 
for the OBS and No Forcing cases (Figure 14). The more absorbing simulations are 
associated with a stronger easterly jet as the dust is leaving Africa (Figures 13 and 17), 
resulting in more rapid transport and generally longer atmospheric residence times (Table 
2). The heating by the dust has a stabilizing effect on the lowest part of the atmosphere 
near the source region but sets up an onshore flow and strengthens the mean meridional 
circulation associated with the Hadley cells (Figure 18). The enhancement of the upward 
motion in the Hadley cells is associated with the peak dust concentrations and leads the 


860 

861 

862 

863 

864 

865 

866 

867 

868 

869 

870 

871 

872 

873 

874 

875 

876 

877 

878 

879 

880 

881 

882 


overall elevation in the dust altitude. Overall we find a better agreement between the 
modeled and observed dust distributions for the simulations with the higher dust 
absorption, leading us to recommend the OPAC or SF optical tables over the OBS tables. 

Finally, the impact of particle shape is a second order effect in our simulations. Although 
particle nonsphericity results in a higher dust mass extinction efficiency and higher dust 
AOT (Table 2 and Figure 2) than for spherical properties, the character of the dust 
transport is more clearly group by refractive index choice (i.e., OBS vs. OPAC vs. SF) than by 
particle shape (Figures 10 and 14). That said, the forcing appears to be sensitive to particle 
shape, with particularly the spheroidal model showing > 10% higher total atmospheric 
forcing over land than the corresponding spherical simulations. Ellipsoidal simulations are 
more similar to their spherical counterparts. We stress that even if radiative effects of 
particle shape are relatively unimportant in the overall dust-climate interaction, as our 
studies here suggest, it is likely an important consideration for observation (i.e., radiance) 
simulation from modeled aerosol fields. 

Future studies with this model would invoke several enhancements to our previous 
modeling capabilities. First, higher spatial resolution model runs would permit a more 
realistic simulation — and, hence, evaluation — of dust from Asian sources. Second, our 
simulations evolved very similar particle size distributions that were insensitive to the dust 
optics and generally underestimated coarse mode dust mass; in future simulations we 
would investigate alternative particle size formulations (e.g., Kok 2011) and a more 
thorough evaluation of the simulated particle size distribution. There is also a need to 
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evaluate model parameterizations, particular with respect to wet loss processes in 
governing the east-west gradient of dust over the Atlantic. Likewise we could as well 
investigate the importance of the wind speed PDF in governing dust emissions, as 
suggested elsewhere (e.g., Cakmur et al., 2004; Su and Toon, 2009). Finally, a coupled 
ocean model was not available in GEOS-5 at the time these simulations were performed; 
such a capability is soon available and can be applied to these same sorts of simulations to 
more fully explore the climate feedback as a function of dust optical properties. 
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Figure 1 . Spectral refractive indices used in this study for the OPAC, SF, and OBS 
simulations. Note that the OBS curve (the orange curve) is identical to the SF curve 
wavelengths greater than 3 pm. 



1152 Figure 2. Spectral aerosol optical thickness (AOT, left), single scattering albedo (SSA, 

1153 middle), and asymmetry parameter (g, right) for each of our optical lookup tables 

1154 considered in this study. Results presented are for the climatological particle size 

1155 distribution at Cape Verde from our OPAC-Spheres simulation normalized to an AOT of 1 at 

1156 550 nm. Absolute values (top) and differences from the OPAC-Spheres values (bottom) are 

1157 presented. 
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Figure 3. 2011 - 2050 climatology of annual mean dust emissions [Tg yr 1 ] for each of our 
simulations. We show the mean (horizontal line), inner quartile range (box), and extrema 
(whiskers) of the time series. 
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Figure 4. Years 2011 - 2050 climatology of seasonal dust emissions for our No Forcing 
simulation. 
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1165 Figure 5. Map of dust surface mass concentration sites summarized in Table 3. Also 

1166 shown is a scatter plot of the OPAC-Spheres modeled climatological mean surface mass 

1167 concentration compared to the multi-year observations at each site. The one-to-one line is 

1168 indicated by a solid line. The dashed lines show the ten-to-one and one-to-ten lines. Colors 

1169 correspond to broad geographical regions to facilitate analysis of scatter plot. 
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Deposition Site Locations 



Figure 6. As in Figure 5, but dust deposition sites listed in Table S2. 
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Figure 7. Years 2011 - 2050 climatology of seasonal 550 nm dust AOT for our No Forcing 
simulation. 
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Figure 8. Climatological JJA aerosol optical thickness from the OPAC-Spheres (shading) 
and No Forcing (contours) model runs (left, dust only) and MISR observations (right). Also 
shown are the locations of the AERONET sites used in this study. 
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1178 Figure 9. Climatological simulated JJA dust AOT at four different longitudes west of 

1179 northern Africa. Also shown is the MODIS Aqua over ocean AOT climatology for JJA 2002 - 

1180 2011. 
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Figure 10. Climatological JJA dust AOT in t 

le latitude band 10° - 30° N for our simulatic 


for the years 2011 - 2050. Also shown is the climatological JJA MODIS Aqua over ocean 
AOT averaged over the period 2002 - 2011. At left the results are plotted on a linear scale. 
At right the model and satellite values are normalized to the same value of AOT at the 
longitude of Cape Verde and the results are plotted on a log-y axis to emphasize differences 
in the simulations. 



1187 Figure 11. Comparison of climatological monthly mean AOT from our simulations to 

1188 multi-year monthly means from AERONET at seven sites impacted by Saharan dust. The 

1189 whisker bars indicate the standard deviation of the monthly mean values. The sites are 

1190 described in Table 4 and their locations shown in Figure 8. The AERONET mean and 

1191 standard deviations are indicated by the grey line and light grey shading. 
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Figure 12. Comparison of JJA climatological simulated dust particle size distribution 
(lines) to the seasonal multi-year AERONET retrieved particle size distribution at four 
AERONET sites. 
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Figure 13. Simulated JJA climatological dust mass mixing ratio profiles at four longitudes 
west of the Saharan source region. Results are shown for three model simulations, as well 
as the year 2011 seasonal attenuated backscatter measured by CALIOP (right column). On 
the model plots we also show the profile of the climatological seasonal mean zonal wind. 
White regions at the surface in the CALIOP plots are below the local topography. 
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Figure 14. Dust median altitude for each of simulations. Results are shown for the JJA 
climatological mean value averaged in the latitude range 10° - 30° N. 
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Figure 15. JJA climatology of TOA shortwave aerosol forcing efficiency. At left we show 
the SF-Spheres and OPAC-Spheres values, as well as their difference. At right we show the 
difference only for the OPAC-Ellipsoids, OPAC-Spheroids, and OBS-Spheres simulations 
relative to the OPAC-Spheres baseline. 
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Figure 16. JJA climatology of difference in mean surface wind speeds (top) and dust 
emissions (bottom) relative to the reference No Forcing simulation. Differences are shown 
for the OPAC-Spheres (left) and OBS-Spheres (right) simulations. 
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Figure 17. JJA climatology of difference in mean zonal temperature profile (shading) and 
zonal wind speed (thin contours) averaged 15° W - 0° against the reference No Forcing 
simulation. The differences are for the OPAC-Spheres simulation (left) and the OBS- 
Spheres simulation (right). The thick contour shows the 100 pg nr 3 dust concentration 
isosurface for the perturbation experiments to indicate maximum dust concentrations. 
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Figure 18. JJA climatology of zonal averaged zonal wind (top) and mean meridional 
circulation (bottom). Contour lines show the values of the reference No Forcing simulation, 
while the shading shows the difference of the OPAC-Spheres (left) and OBS-Spheres (right) 
relative to the reference simulation. The thick contour shows the 20 pg nr 3 dust 
concentration isosurface to indicate maximum dust. The arrows give a qualitative sense of 
the vertical transport. For the mean meridional circulation, negative values indicate 
counterclockwise transport (contours) or an increase in the counterclockwise transport 
(shading). 


